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l/~j Abstract 

^^ We study mathematically a system of partial differential equations arising in the 

n , modelling of an aging fluid, a particular class of non Newtonian fluids. We prove 

.^rfH well-posedness of the equations in appropriate functional spaces and investigate the 

. longtime behaviour of the solutions. 
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1 Introduction 



Our purpose is to study mathematically a system of partial differential equations arising 
in the modelling of some particular non Newtonian fluids. These fluids are often called 
00 o-ging fluids. Two physical phenomena are indeed permanently competing within the flow 

^^ of such fluids. On the one hand, the fluid ages in the sense that it solidifies. On the other 

Q hand, aging is counterbalanced by a flow-induced rejuvenation. 

^J The specific modelling we consider has been proposed in j3| on the basis of phenomeno- 

(^ logical arguments and experimental observations. A coefficient /, called the fluidity en- 

CN codes aging for all times and at every location within the fiuid. The fluid is solid where 

/ = 0, and behaves all the more as a liquid when / grows. Our mathematical study aims to 
. ^ contribute to better understand how well such a model captures the essential phenomena 

^ at play in fluid aging. 

For our study, we proceed in a one-dimensional setting corresponding physically to the 
consideration of a laminar Couette flow. Our three unknown fields, the velocity u, the 
shear stress r and the fluidity / are deflned as functions of a space variable y varying 
in the interval [0, 1]. They are also, of course, functions of the time i > 0. The specific 
system we choose for our study reads 

du d^u dr 



dl 
dt 



{-l + ^\T\)f-uf. (l.lc) 



Six dimensionless coefficients, all positive, constant in time and throughout the domain, 
are present in the system: the density p, the viscosity r/, the characteristic relaxation 
time A, the elastic modulus G, and two coefficients ^ and v specifically related to the 



equation for the evolution of the fluidity /. System (1.1) is a fully coupled system of 
three equations. The flrst two equations are classical in nature. The first one is the 
equation of conservation of momentum for u. The second equation rules the evolution of 
the shear stress r. The non classical ingredient therein (as opposed, say, to an Oldroyd-B 
type equation) is the presence of an extra parameter, the fiuidity /, the role of which 
is formally similar to that of an inverse time in a relaxation phenomenon. We note that 
when / is a constant in time and throughout the domain, the equation agrees with the one- 
dimensional Oldroyd-B equation considered e.g. in [21 [5]. The third equation is of the form 
of one of the many such evolution equations suggested in [1]. We hope it is, in this respect, 
prototypical of a general class. It models the evolution of the fiuidity / in function of the 



stress tensor. The right-hand side of (1.1c) may differ from one model to another. The 
important ingredient is the presence of two competitive terms: a negative term modelling 
aging and a positive term modelling rejuvenation. For mathematical convenience, we have 
taken two particular instances of these two terms. 



We examine well-posedness and longtime behaviour for system (1.1). Because we 
provide self-contained proofs, our study is rather long. Similar questions on a different, 
although related model for a viscoelastic fiuid, have been examined in [QlfTO]. 

Our article is articulated as follows. 

To start with, we prove in Section 2 that the system under consideration admits a 
global-in-time solution in appropriate functional spaces. The solution is shown to be 



unique, and indeed strong. System (1.1) is thus satisfied in a classical sense. Our precise 
statement is the object of Theorem 2.1 The bulk of Section 2 consists of our proof. The 
arguments are standard arguments of mathematical fluid dynamics: formal a priori esti- 
mates, approximation, rigorous a priori estimates, convergence. The many nonlinearities 



present in system ( 1.1 ) however prevent us, in the current state of our understanding, from 
extending our analysis to settings in dimensions higher than or equal to two. Technically, 
this is related to the fact we repeatedly use, in our arguments, that H^ functions are L°° 
functions, a specificity of the one-dimensional setting of course. 

In Sections 3 and 4, we study the long time behaviour of the solution. Section 3 deals 
with return to equilibrium. We supply the system with homogeneous Dirichlet boundary 
conditions for the velocity and investigate whether the flow converges to a steady state. For 
homogeneous Dirichlet boundary conditions, the steady states are (u = 0,r = c, / = 0), 
where c is a constant throughout the domain. The long-time convergence to these steady- 



states sensitively depends, in system (1.1), of the fluidity /. The situation is qualitatively 
different depending on the fluidity /o at initial time. The more delicate, but of course 
more interesting, case mathematically is the case where the fiuidity /o at initial time does 
not vanish everywhere: a part of the material, possibly the whole of it, is originally fiuid. 



Section 3.1 addresses this case. We show (and the proof is quite substantial even in the 



one-dimensional setting we consider) that the fiow converges to the null steady-state in 



suitable functional norms. The precise statement is the purpose of Theorem 3.1 The 



convergence is then shown to be polynomial in time, for all three fields u, r and /. The 



rates of convergence are made precise in Theorem 3.2 Numerical simulations we perform in 



Section 5 will show these rates are indeed sharp. It is interesting to emphasize the physical 
signification of our mathematical results. With regard to modelling, the convergence of 
the fluidity / to zero that we establish, under homogeneous Dirichlet boundary conditions, 
means that when left at rest, the fluid progressively solidifies, a certainly intuitive fact. In 
addition, for u and r, the rate of convergence sensitively depends on the size of the region 
where, originally, the material is liquid (a size measured by our parameter /3 defined in 



(3.17) and present in the right-hand sides of the estimates of Theorem 3.2). The larger 
the liquid region the quicker the convergence of both the velocity and the shear stress to 
zero. It is not completely clear to us whether the latter qualitative behaviour is or not 
compatible with experimental observations or physical intuition. 

If the material is entirely solid at initial time, that is fo = everywhere, the behaviour 
is quite different. Then the material stays solid for all times, while the velocity and 
shear stress vanish exponentially fast. We present the simple analysis of this behaviour in 



Section 3.2 Note that the result agrees with simple physical intuition. 

Non-homogeneous boundary conditions, studied throughout Section |4j are, as always 
for questions related to long-time behaviours, significantly more intricate to address. We 
adopt constant boundary conditions, respectively u = and u = a>Oaty = and y = 1. 



We begin by showing in Section 4.1 that, when we impose that the fluidity is strictly 



positive everywhere, there exists a unique steady state. We next show in Section 4.2 
that this steady state is stable under small perturbations, our precise result is stated 
in Theorem |4.1[ When the perturbations of the state are not small, analyzing return to 
equilibrium is, in general, beyond our reach. We are however able to show that, when we 
assume a particular form of the initial condition (namely linear velocity, constant shear 
stress, constant positive fluidity), then return to equilibrium does hold true even if the 
initial condition is not close to the steady state. Some suitable assumptions relating the 



size of the parameters in system (1.1) and the non-zero boundary condition a are also 



needed (see condition (4.21 )). Our precise result is Theorem 4.2 The reason why we have 
to assume this specific form of the initial conditions is purely technical (and our numerical 
simulations will actually show that these restrictions are, in practice, unnecessary). In that 



case, system (1.1) reduces to a two-dimensional system of ordinary differential equations, 
for which Poincare-Bendixson Theory allows us to understand the longtime behaviour. 
Our study is performed in Section |4.3[ 

As briefiy mentioned above. Section 5 presents some numerical simulations. We first 
show that the rates of convergence estimated by our various mathematical arguments in 
the various regimes considered in Sections [3] and |4] are indeed sharp. We also investigate 
numerically the stability of the steady state. Our simulations show that, irrespective of 
the size of the initial perturbation (and thus in a more general regime than that for our 
mathematical arguments), the fluid returns to equilibrium, or more generally converges to 
the suitable steady state. The rates of convergence are also examined. 

We conclude this introduction by mentioning that, despite their limitations, our results 
show that the model derived in j3] does adequately account for aging and rejuvenation. 
However, two shortcomings need to be emphasized. Both originate from the mathematical 



nature of equation (1.1c) (and are actually related to the fact that the Cauchy Lipschitz 
theory applies to this equation). First, when / vanishes, then / remains zero for all 
subsequent times. This property, present everywhere in our mathematical study, pre- 



vents fluidification to occur after solidification. Tliis clearly limits the range of materials 
covered by the modelling (compare muds and concrete, say). Second, / can only van- 
ish asymptotically and never in finite time unless it is already zero before. Otherwise 
stated, solidification can occur, but never in finite time: again a modelling limitation. 
The one usefulness, if any, of our study, is therefore to point out that a mathematically 
well founded model where fluidification and solidification compete on an equal footing is 
still to be derived. Our study implicitly points to suitable directions to this end. 

Further mathematical investigations on models for aging fluids will be presented in [1] . 



2 Global existence and uniqueness 

In this section, we establish the following global existence and uniqueness result for sys- 
tem ( |1.1[ ) supplied with initial conditions uo,to, /o and the boundary conditions u{t, 0) = 
and n(t, 1) = a > for all time t € [0, T] (where a is a constant scalar). 

Theorem 2.1 Recall that J7 is the one-dimensional domain [0, 1] and that T > is fixed. 
Consider the initial data 



{uo,ToJo)eH^{nf withfo>0. 



(2.1) 



Then there exists a unique global solution (n, r, /) to system ( 1.1 ) such that for any T > 0, 

iu,T,f)e{Ci[0,T];H')nL\[0,T];H^))xCi[0,T];H^)xC{[0,T];H') (2.2) 

and / > for all x & 0, and t G [0, T] . 
In addition, we have 

' du dr df ' 



€L' 



,r];L2)xC([0,r];L2)xC([0,r];L2), 



(2.3) 



so that the equations in (1.1) are all satisfied in the strong sense in time. 



Before we get to the proof, we eliminate the non-homogeneous Dirichlet boundary 
condition, introducing the auxiliary velocity field 

u{t,y) -ay. 

This velocity field, which we still denote by u, solves the system 



du 
~di 
dr 

m 

dl 
dt 



d^ 



dr 



dy"^ dy ' 



G 



du 



fr + Ga, 



dy 

(-1 + CM)/' 



^f, 



(2.4a) 
(2.4b) 

(2.4c) 



supplied with the homogeneous Dirichlet boundary conditions u{t, 0) = and u{t, 1) = 
for all time t £ [0, T] and initial conditions uo,to,/o S H^{Q). The proof of Theorem 



will actually be completed on system (2.4). The result on (1.1) then immediately follows 



2.1 



Proof. The proof falls in eight steps. The first five steps consist in deriving formal 
a priori estimates. These estimates are next made rigorous for a sequence of approxi- 
mate solution in Step [6j The convergence of this sequence is proven in Step [7j thereby 
establishing existence of a solution to (2.4). Step [S] addresses uniqueness. 



Step 1: Non-negativity of the fluidity. Let us first formally prove that / > 0. 
Fix y G il and introduce 

Eo = {y £ njoiy) > 0} . 

For y £ Q\ Eq, we have fo{y) = and thus f{t,y) = for all time t £ [0,T] because 



of (2.4c). On the other hand, for y G Eq, we now show that f{t,y) > for all time t G 



[0,T]. We argue by contradiction and suppose, by continuity of f{-,y), that 

i^ = inf{tG(0,r],/(i,y)=0}<r. 



The Cauchy-Lipschitz Theorem applied to ( 2.4c| ) with zero as initial condition at time tm 
implies that /(t, y) = for t G (tm — ^,'tm + e) for e > 0, which contradicts the definition 

of tm- 

We have therefore shown that / stays zero where it is zero, and stays positive where 
it is positive, which in particular implies non-negativity everywhere. 

Step 2: Formal first energy estimates. We again argue formally. We first mul- 



tiply the evolution equation (2.4a) on u by u itself and integrate over the domain. This 
gives a first estimate 






du 
dy 



(t, 



L2 Jn \dy 



(2.5) 



Similarly, we multiply the evolution equation (2.4b) by r and integrate over Q, to find 

' du 



^A|lk(V)lli.+ 



Mt, 



.' Gy^^-.)(,.) + G..(t). 



where we denote by 



q{t) = I q{t,y)dy 
In 



(2.6) 



(2.7) 



the average over fi of a function q : (t, y) G [0, T] x 0, 



Combining estimates (2.5) and (2.6) and using integration by parts and the fact that u 

du 



vanishes on the boundary, we obtain 

[Gp\\u{tr)\\h + X\\r{t,-)\\h)+Gv 



1 d_ 
2dt 



dy 



(t, 



+ 



L2 



^/fT){tr) ^=Gaf{t). 



L2 



We now turn to (2.4c). Integrating (2.4c) over Q. yields 



(2. 



(2.^ 



The Young inequality 



^r\f = ^f^ 



\T\P<^-f + ^fT^ 



,1 , i^ 



2v' 



then yields 



dt 



ll/(*,-)|lLi + ||/(t,-)lli. + ^ll/(t,-)lli3<|^||(y7r)(t,-) 



2 

L2 



(2.10) 



Collecting (2.8) and (2.10), we obtain 



2dt 



(^GpMt,-)\\h + X\\T{t,-)\\l, + '^^\\f{t,-)\\j 



+ Gri 



du , 
dy 



it,' 



+ 



L2 



(y7r)(t,-)|[, <C7a||T(t,.)||i2 



(2.11) 



where C, here and throughout our text, denotes a constant, the actual value of which is 
independent from T and only depends on the domain Q and the coefficients p, rj, A, G, ^, u 



in (2.4) . 



Applying the Gronwall Lemma to (2.11), we obtain 

sup (\\u{t,.)\\l, + \\r{t,.)\\l, + \\f{t,.)\\^A+ r (\\n{t,.)\\l,+ (/7r)(t,.) 
<e[o,T] ^ '' Jo ^ ^ ^ 

< Cot, 



L2 



dt 



(2.12) 



where Co,t is a constant depending not only on Q, p, -q, A, G, ^, z^, but also on the boundary 
condition a, the initial data uq, tq, /o and the time T. 



Remark 2.1 For homogeneous boundary conditions, that is a = 0, we mention that the 



right-hand sides of (2.8) and thus (2.11) vanish. The constant Cq^t in (2.12) therefore 
does not depend on T and we get a bound uniform in time. 



Step 3: A priori estimates on an auxiliary function. Denote 



g{t,y)= / {r{t,x)-fit))dx. 
Jo 

This function g satisfies homogeneous Dirichlet boundary conditions and formally solves 

d'^g dr 
dy"^ dy 



Using (2.4a) and (2.4b), which respectively imply 



du d'^ 



Pi^ = Vt^ {u+z9] , 



1 

— i 



and 



dg 



f\fT-fT)dx + Gu, 

Jo 



we remark that the auxihary function 



solves: 



1 fy 

U = U+ - (T-f) 

V Jo 

1 

= u+ -g 
V 



dt p (92/2 Ar/ 7o ^^ 



(2.13) 



(2.14) 



Multiplying equation (2.14) by 






and integrating over yields 



ld_ 
2dt 



dy 


L2 P 


dy^ 


%■) 


2 

L2 


<c(||(/ 


r){t,.)\\ 


-I 

Jn 


dy^ 


( 



Jn 



dy"^ 



it 



' / / ' 



after elementary manipulations in the right-hand side. Then using the Young and the 
Cauchy- Schwartz inequalities, we obtain 



1 d^ 

2dt 



dU 

dy 



it,' 



L2 2p 






2 
L2 



<c(\\f{t,■)\m{^/fr){t, 



) ^, + Mt,-)\\h 



(2.15) 



Using (2.12), we know that the right-hand side is L^(0,T). In view of our regularity 



assumptions on the initial conditions, we have 



dU 
dy 



therefore deduce from (2.15) that 



,.o = lr^J<^°-^'^^'<">'*' 



U e L°°i[0,T],H^)nL'^{[0,T],H^). (2.16) 

Step 4: L°° estimates. We are now in position to obtain (again formal) L°°-bounds 



on r and /. We consider the evolution equation (2.4b), which we rewrite in terms of U 



defined by (2.13) and using r defined for r as in (2.7): 



dt dy 



T] J ri 



Multiplying this equation by r, we obtain 



A d_ 
2di' 



/ + 



G 



kr < C Irl- 



dU 



dy 



+ rl ■ Irllia +a|r 



so that, repeatedly applying the Young inequality, 



A d, ,2 f, G 



kr < c 



dU 



dy 



I M-< I 



(2.17) 



dU 



We apply the Gronwall Lemma to (|2.17D and use ^^ G L^([0,T],L°^) because of (2.16) 

dy ^ ' 



estimate (2.12) and tq G H^{Q) to obtain 



Tt.) 



IL°° 



<Co' 



(2.18) 



-f-'yf)f + ^\r\f, 



that is, r GL°°([0,r],L°°). 

As for the function /, using the Duhamel formula for the evolution equation (2.4c) 
rewritten as 

dt 
we obtain, for almost all y £ Q, 

f{t,y) = e-/o(/+^/^)(^>^)'^7o(y) + e re-r(/+^/^)(*''^)'^*' {\r\f) {s,y)ds 

Jo 

< h{y) + - WrU^iL^, f e-^>f'^''^y>^''vf\s,y)ds, 



V -tI^-);, 



where we have used the non-negativity of / and the previously derived L°°-bound on r to 
obtain the second line. The above equation leads to 



/(t,y) < /o(y) + \ l|r|Loo(^^) (i - e-/o -/^(«.^)'^^) 



< /o(y) + -, l|r|L, 



y -^riL^) 



(2.19) 



Using that /o G H^ and that we work in a one- dimensional setting, we obtain that / G 

L°°([0,r],L°°). 



Remark 2.2 For homogeneous boundary conditions, the Gronwall Lemma applied to (2.17) 
implies 



|r(i,-)|lioo<||ro||i2e ^v' + 



dU 
dy 



(s,-) 



ds + sup ||T(t, •)||^2 / e^-i 
L°° te[o,T] Jo 



(s-t) 



ds. 



(2.20) 



Moreover, as explained at the end of Step^ the constant Cq^t in (2.12) does not de- 
pend on T. Hence, the right-hand side of (2.15) and the bound in L'^([0,T]; L°°)-norm. 
dU 



for ——, deduced from (2.16), also do not depend on T. It follows from (2.20) that the L 
dy ^ ' 



bound (2.18) on t is uniform in time. Equation (2.19) yields a similar conclusion for the 



bound on f . 



Step 5: Second a priori estimates. In order to get estimates on higher order 



derivatives, we now differentiate with respect to y the evolution equation (2.4b) and obtain 



d fdT\ _ ^d'^u dr df 
dt \dy ) dy"^ dy dy 

- c— - -— _ f^ _ ^ 
dy"^ T] dy dy dy 



(2.21) 



Likewise, we differentiate witli respect to y the evolution equation (2.4c) and get 
dt \dy ) dy dy dy ' 



(2.22) 



dr 



df 



Multiplying equations (2.21) and (2.22) respectively by — - and — — , integrating over the 

oy oy 



domain, summing up and using that both r and / are in L°°([0,T],L°°), we obtain 



1 d 

2dt 



dr 
dy 



(t,-) 



<Co' 



+ 
d^Udr 



dl 
dy 



it,- 



L^ . 



9r\ df dr d\T\df fdf\ 



dy"^ dy \dy J dy dy dy dy \dy J 



(t,- 



Repeatedly applying the Young inequality and using that r belongs to L'^{[0,T],H^] 



(which implies 



d\r\ 



dy 

1 d j 
2dt I 



{t,y) 

dr 



dr 
dy 



it,y) 



X 



dy 



it,-) 



+ 



L2 



<C7o- 



dr 
dy 



it,- 



for almost all t £ [0,T] and y G O), we obtain 
d^U 



dy 



+ 



L2 



dl 
dy 



it, 



+ 



L2 



dy- 



-it,' 



2 ^ 



2.23) 



We apply the Gronwall Lemma to (2.23), use that tq, /o G H^{Q,) and the estimate (2.16) 
to obtain that r, / G L°°([0,T], i^^). 



It follows from the definition ( |2.13D and the estimate ( |2.16D that u G L°°{[0,T], H'^) n 
L\{0,T],H^). 

Step 6: Construction of an approximate solution. Now that we have estab- 
lished all the necessary formal a priori estimates, we turn to the construction of a sequence 



of approximating solutions to (2.4) on which we will rigorously derive these a priori esti- 



mates. We introduce, for n > 1, the sequence of systems 

dur,. 



P- 



dt 



V- 



d'^Un dT„. 



dy"^ dy ' 



dTn dUn 
A^T- = G^- Jn-lTn + txO, 

dt dy 



dfn 

, dt 



i-l + C\Tn\)fn-lfn-iyfn^lfl 



(2.24a) 
(2.24b) 
(2.24c) 



supplied with the homogeneous Dirichlet boundary conditions u„(t, 0) = and Unit, 1) = 
for all time t G [0,r] and initial conditions {uno,Tno, fno) = iuo,To, fo)- We actually use 
the initial condition (uoi''"Oi/o) also to initialize the iterations in n, thus the coincidence 
of notation. 

We argue by induction. Consider 

(u„_i,r„_i,/„_i) G iC{[0,T];H^)nL\[0,T];H^)) x C{[0,T];H^) x C{[0,T];H^) 



and fn-i > 0. We first show that there exists a unique solution {un,Tn,fn) to (2.24) 
belonging to the same functional spaces and such that fn > 0. For this purpose, we 



decompose (2.24) into two subsystems: the linear (Oldroyd-B) type model coupling the 



evolution equations (2.24a) on n„ and ( 2.24b[ ) on t„ on the one hand and the ordinary dif- 
ferential equation (2.24c) on fn satisfied for ally G il. on the other hand. The existence and 



uniqueness of a solution {un, t„) in the space (C([0, T];H^) D L'^{[0, T];H^)) xC([0, T];H'^) 
for the former system is obtained using a classical approach (see for instance [^ for a very 
close system). We now turn to /„. We show that /„ exists in C{[0,T];H^) and /„ > 0. 
The equation ( |2.24c ) writes 



dfn 



dt 

Jn\t=0 = /o 



i^{t,fn,y), 



(2.25) 



where ^ is a function from [0, T] x M xil to M. 

We first fix y S rj and show that the function fn{-,y) is continuous in time and non- 
negative. The function ip is continuous in its first two variables and locally Lipschitz 
in its second variable. The Cauchy-Lipschitz Theorem shows there exists a unique local 
solution with fo{y) as initial condition. Let [0, T*) be the interval of existence of the 
maximal solution for positive time. For all t G [0,T*), we have /„ > 0, using Stepfl] In 
addition, since fn-i and /„ are both non-negative, ( |2.24c ) implies for all t € [0,r*), 

dfn 

dt 



— S rn\ Jn—ljn 



< ^ llr I 
_ ? II 'ra| 



Ct{L°°) \\jn-l\\cT{L^) 



Jni 



(2.26) 



using that both t„ and fn-i belong to C{[0,T];H^). The Gronwall Lemma then proves 
that fn remains bounded on [0, T*] and thus we have established existence and uniqueness 

on [o,r]. 

We now turn to the local property of continuity of /„ as a function of y. We use that 
the function ip is continuous in y, because both t^ and fn-i are continuous in y in our 
one-dimensional setting and the theorem on the continuous dependence on a parameter 
for ordinary differential equations of the form ( 2.25| ) (see e.g. [31 Theorem 1.11.1, p. 126]). 

^ hp1nTio-« tn n(\r\ T]-T'^'^ 
dy 

dfn 

following linear ordinary differential equation on — — 

dy 



We now show that — — belongs to C([0,T];L^). We consider, for almost all y G J7, the 



d^fdfn 
dt \ dy 

where we have introduced the functions 



A^-^ + B, 
dy 



A = i\Tn\fn-l 



fn-l — 2z^/«_l/n 



B = i- 



dy 



f f \ 1 1\ I ^\ t djn-1 f2djn~l 

fn-lfn + {C\Tn\ - 1) fn^ I'fn 



dy 



dy 



eC([o,r];L°°), 

gC{[0,T];L^). 

dfn 



(2.27) 

(2.28) 
(2.29) 



The Cauchy-Lipschitz Theorem then guarantees the existence of — — (•,?/) continuous in 



time, for almost all y G fi. The Duhamel formula applied to (2.27) yields, for all t G [0,T] 



10 



and almost all y G fi, 



dfn 



dfo 



dfn 



so that, using (|2.lD, i\2.28\) and (\2.29\i, ^ belong to C([0,r];L2). As ( |2.27D is the 



derivative with respect to y of (2.24c), this yields /„ G C{[0,T]; H^). 



Now that we have established, for all n, the existence of a solution (u„, t„, /„) to (2.24) 



in the appropriate functional spaces (as in (2.2)-(2.3)), we derive, for {un,Tn, fn), the a 



priori estimates formally established on {u,t, f) in the previous steps. Estimate (2.8) now 
reads 



1 d 
2dt 



[Gp\\un{tr)\\h + X\\r4t,-)\\l,) + Gr^ 



dUn 



dy 



(t,-) 



+ 



L2 



(\//n-lT„j (t, •) 



2 
L2 



Ga7v;:(t). (2.30) 



Likewise, (2.10) is now replaced by 



d^ 

dt 



\\fn{t,-)\\L^+ (/„-l/„)(t,-) + ^/ ifn^lf^) it,-) < ^\\{y^U^irn) (t,-) 



2 
L2 



(2.31) 



Collecting (2.30) and (2.31) yields the following estimate, analogous to (2.11), 
~ (^Gp ||n„(t, ■)\\l, + A ||r„(t, ■)\\l, + ^ ||/„(t, OUz^i 



+ Gr] 



dUn 



dy 



it,- 



+ 



1 



L2 



[^/f^lrn) (t, oil', < Ca ||r„(t, OH^a , (2.32) 



and therefore, (2.12) holds with {un,Tn,fn) instead of {u,T,f). 

The arguments given in Step [3] to derive (2.15) and in Step [4] for the L°° estimates can 
be mimicked for the approximate system in {un,Tn, fn-i) instead of {u,T,f), and the 
corresponding auxiliary functions gn and C/„. 

At this point, we have rigorously established on {un,Tn, fn) and our formal estimates 
of steps [2] to [4) 



sup sup (||n„(t, •)|Il2 + ||r„(t,-)|lL2 + ||/n(i, OIIlO < Co,r, 
" te[o,T] 



(2.33) 



and 



sup sup {\\Un{t, ■)\\hi + \\Tn{t, OIL- + \\fn{t, OIL-) + l|f^nlL2 (^2) < Co,T, (2.34) 

n te[0,T] -' 

where we recall that Co,t denotes various constants which depend on the coefficients in 



system (2.4), the initial data uo,to,/o and the time T. 
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We now turn to the a priori estimates of Step [5J Using arguments similar to the formal 
arguments of Step [5j we obtain 



dt 



dxn 
dy 



<Co, 



T 



(t,-) 



dy 



+ 



L2 



it, 



dfn 

dy 



it,' 



L'2. 



+ 



L2 



dfn- 



dy 



Ht, 



+ 



L2 



dy 



-it, 



+ 



L2 



d^Un 



dy^ 



it,-) 



(2.35) 



^2, 



We now observe that showing H^ bounds on r„ and /„ is less straightforward than in our 
formal Step [5| We introduce 



Ynit) 



dTn 

dy 



t,-] 



+ 



L2 



dfn 

dy 



it,-) 



L2 



Integrating (2.35) from to i < T, we see that Yn satisfies 



Ynit) < Co,T I Yn + Co,T [ Y„,-l + Co,T HC/nH^ ,^2^ + Yq. (2.36) 

Jo Jo ^ 



Applying the Gronwall Lemma to (2.36), we find 



Ynit) < (Co,T||C/n|li|(j^2)+yo) 6^°.^* + Co,T / l^-l (s)e^«'^(*-^)d6 

I' u 

which we rewrite 

Ynit) < Co,T + Co,T / Yn-lis)ds. 

Jo 

Arguing by induction, one can check that this implies, for all t G [0, T] and n. 



n-l 



iyO.c.E^ + ^^^^'o. 



i=0 



n! 



It follows that, Co,T denoting various constants. 



sup sup Ynit) < Co,re^«'^^, 
n te[o,r] 



(2.37) 



Recalling that ^i^ — = —;^ (t„ — r^), we use inequalities (2.34) and (2.37) to derive 



dy dy i] 

sup sup i\\Unis,-)\\Hi + \\Tnis,-)\\Hi + \\fnis,-)\\Hi) + \\Un\\L2fH2)<Co,T- (2.38) 



n se[0,T] 

This implies 



sup sup 

n se[o,r] 



dUr. 



dt 



is,-) 



+ 



LliL^) 



dTn 

dt 



■is,-) 



+ 



L|(L2) 



dfn 

dt 



is, -) 



< Co,T- (2.39) 
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Step 7: Convergence of the sequence of approximate solution. The bounds 



of a subsequence, we have the weak convergences 



obtained in the previous steps, namely ( 2.38 ) and ( 2.39 ) show that, at least up to extraction 
lave the w( 

yUm Tm In 
'dUn dTn dfn 



(n, r, /) weakly- * in L~([0, T]; H 



1^3 



Ur, 



u weakly in L 



,T]:H' 



(2.40) 
(2.41) 



dt'dt' dt 



fdu dr df\ . 2, 



,T]:L 



2\3 



(2.42) 



But, in order to pass to the limit in (2.24), we need the convergence of the whole sequence 



itself because (2.24) involves indices n — 1 and n and strong convergence to establish 
convergence of the product terms fn-ifn, \Tn\fn-ifn, fn-ifn- 

We now establish strong convergence of the whole sequence. We prove this convergence 
in (-L~([0, T]; L'^{n)))^. This wiU a posteriori imply that ah the convergences ( |2.40[ ), ( |2.41[ ) 
and ( |2.42 ) actually hold true not only for a subsequence, but the whole sequence itself. 
And this will provide sufficient information to pass to the limit in our nonlinear terms. 

We introduce the notation: hn = hn — hn-i and derive the evolution equations for 

\Um Tni Jn) 



dUn d'^Un , dTn 



dt 



7] 



dy^ dy 



dt dy 



Jn—lTn ''"n— l/n— 1) 



dfn 

dt 



-1 +^|t„_i|)(/„_i/„ + /„_i/„_i) 

- ^fn-l{fn + fn-l)fn - '^fn-lfn-1 + C\Tn\fn-lfn- 



(2.43a) 
(2.43b) 

(2.43c) 



Since {un,Tn,fn) belong to the spaces that appear in (2.2) and (2.3), the same holds 



for {un,Tn, fn)- We multiply equations (2.43a), (2.43b) and (2.43c), respectively by Un,T, 



and fn, integrate over Q, sum up and use the non-negativity of fn-i and fn to find 



±{Gp\\Unit,-)\\l2+X\\Tn{t,-)\\l2 + 



fn{t, 



L2 



< 



Tn-lfn-lTnit,-) 



n 



+ / (-l + Ckn-l|)(/n-l/„ + fn-lfn-lfn) - J^fn^lfn-lfn + (,fn-lfn\rn\fn{t, ■)■ 

In 
The presence of two indices n — \ and n again makes an additional step necessary. We 

2 



introduce Xn{t) 



I^n(i,-)|li2 + ||r„(t,-)||i2+ fn{t,-) 



L2 



. Repeatedly using the L° 



bounds (2.34) on {t„, /„,r„_i, /„_i} and the Young inequality, we see that Xn satisfies 

Xn{t) < Co,T(^n(t) + Xn-l{t)). (2.44) 



Applying the Gronwall Lemma to ( 2.44[ ), we find 

Xn{t) < Co,T [ X„_i(s)e^o-^(*-^)(is < Co,Te^°^^^ / Xn^i{s)ds, 
Jn Jq 
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which imphes that 



Xn{t) < — sup Xi{s). 



(n-1)! 



se[o,T] 



^T];L^{n))f. The 



The sequence {un,Tn,fn) is therefore a Cauchy sequence in (L"* 
sequence converges in this space. 

Now that we have strong convergence of the whole sequence, we show how to pass to the 



hniit in all the terms of (2.24), including the nonlinear ones. We only consider \Tn\fn-ifn- 



The other terms can be treated using similar arguments. We use a classical compactness 



result |8i Theorem 5.1, p. 58] to deduce from (2.40) and (2.42) that r^ and /« strongly 
converge respectively to r and / in L^([0, T]; L^)"^. Moreover, fn-i strongly converges to / 
in L°^([0,T];L2). We thus have convergence for |r„|/„_i/„ in L^{[0,T];L'^). 



The triple {u,T,f) thus satisfies system (2.24), at least in the weak sense. We now 
derive further regularity. We have 

uGL\[0,T];H^) with ^ ^ l\[0,T];L^), 



and therefore, by interpolation (see [121 Chapter 3, Lemma 1.2]), 

ueC{[0,T];H^)nL\[0,T];H^). 



Moreover, we have 



dr 9/ 
at' 5t 



2\2 



GL'i[0,T];L') 



and, using the second a priori estimate (2.23) , 



d dr d df 



2\2 



GL'i[0,T];L') 



so that, 



dt dy ' dt dy ^ 

(r,/)GC([0,r];Fi)2. 



We have obtained (2.2) and therefore (2.3), using system (2.4). The non-negativity of the 



fluidity is preserved, passing to the limit. This completes the existence proof. 



Step 8: Uniqueness. Consider (ni,Ti,/i) and (n2,T2,/2) satisfying (2.2 



and so- 
lutions to system (2.4) supplied with the same initial condition {uo,To,fo) G H^{Cl). We 
introduce {u = U2 — ui,t = T2 — Ti, f = f2 — fi) which therefore satisfies 



(2.45a) 
tJ, (2.45b) 

(2.45c) 



du 


d^u 


Bt 


'dt- 


-"V 


dy 


Bt 


„du 




A^ = 


-G ^ 


- Ut 


dt 


dy 




df 
<dt 


-(/i + 


/2)/ 



f2)f-l^{f! + flf2 + fl)f, 
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supplied with homogeneous boundary conditions and (0, 0, 0) as initial data. Multiplying 



equations (2.45a), (2.45b) and (2.45c), respectively by u, r and /, integrating over $7, 



summing up, using the L°°-bounds established in Step [4] for terms involving ti,T2, fi, f2 
and repeatedly applying the Young inequality, we find 



1 d_ 
2dt 



pG\\u{t,-)\\l2 +X\\T{t, ■)\\l2+ fit,-) ^, <Co,T ||r(t. 



\h + 



fit, 



The Gronwall Lemma then implies uniqueness. This concludes the proof of Theorem 2.1 







Longtime behaviour for homogeneous boundary condi- 
tions 



In this section, we study the longtime behaviour of system ( 1.1 ) supplied with homogeneous 



boundary conditions. We will show convergence to a steady state and establish a rate for 
this convergence. For homogeneous boundary conditions, the //^-steady states of (1.1) 
such that / > are exactly the states (u = 0,t = c, / = 0), where c is a constant 
throughout the domain. 

Indeed, such a steady state (uoo, Too, /oo) satisfies, combining equation (1.1a) integrated 



over the domain and (1.1b), 



/oo + l 



(3.1) 



where c is a constant over the domain. We now distinguish between two cases. Ei- 
ther c = 0, in which case Too = 0. The homogeneous boundary conditions on u and (1.1b) 
imply that Uoo = and (1.1c) that /oo =0. Or c 7^ and it follows from (3.1) that Too 

dUoo 



is non-zero and has a constant sign and from (1.1b) that 



dy 



has a constant sign. Be- 



cause of the homogeneous boundary conditions on the velocity, we obtain that Uoo 
Therefore, (1.1b) yields /oo^oo = and /oo = 0, because Too is non-zero in this case. 



0. 



We will show that the longtime behaviour differs both in terms of steady state and 
rate of convergence, depending whether /o ^0 or /o = 0. When /q ^0, a case studied in 



subsection 3.1 , the solution (n, r, /) converges to the steady state (0, 0, 0) in the longtime 
and the rates of convergence are power-laws of the time. In the case /o = 0, the fluidity / 



vanishes for all time, as easily seen on ( 1.1c). In subsection 3.2, we show that [u, r, /) then 
converges to (0,7^,0) in the longtime at an exponential rate, 7^ being the average of tq 
over 0. Evidently, the former case /o ^0 require more efforts than the latter case /o = 
where / = for all times. 

3.1 Case/o^O 

In this subsection, we consider the case /o > 0, /o ^0. We first establish the convergence 
in the longtime. 



15 



Theorem 3.1 Supply system ( 1.1 ) with homogeneous boundary conditions and initial con- 
ditions that satisfy (2.1) and /q ^0. The solution {u,T,f), the existence and unique- 



ness of which have been established in Theorem 2.1, converges to the steady state (0,0,0) 
in H\Q}) X L°°(0) x L°°(0) in the longtime: 



\<t,-)\\H, + \\T{t,.)\\^^ + \\f{t,.)\\, 



0. 



Proof. The proof falls in three steps. In the first step, we establish a lower bound 
for the average of the fluidity /, which, in Step [2j is useful to prove convergence in 
the longtime in L^(Q). In the third step, we show convergence of (^(t, •), T(t, •), /(t, •)) 
in H^{^) X L°°(fl) X L~(0). 

In this section, Cq denotes various constants that are independent from time, while Cj, 
i=l,...,4 denote some fixed constants independent from time. These constants Cq and Ci 
used to be denoted Co,t in the previous section. The subscript T is now omitted because, as 



explained in Remarks 2.1 and 2.2, the constants are independent from T for homogeneous 
boundary conditions. 



Step 1: A lower bound for the average of /. We first derive a lower bound 
on /, defined as in (2.7), and not directly on / because the latter may vanish (since /o 
may vanish) on some part of the domain. Since /o ^0, there exists, by continuity of /o 
(assumed in H^), a non-empty closed interval ilo in $7 where /o does not vanish. Arguing 
as in Step [T] of the proof of Theorem 2.1, the fluidity / does not vanish for all t > 
and y G Oq. The evolution equation (1.1c) on / rewrites, for all t > and y G VLq, 



Wt~f 



l-^\r\+vf. 



(3.2) 



As explained in Remark 2.2, the L°°-bounds on r and / are uniform in time for homoge- 



neous boundary conditions. The equation (3.2) thus implies, for all y G f^o and f > 0, 

d 1 



dtf 



<Co, 



and therefore. 



f{t,y)> 



> 



1 

1 



L°°(Oo) 



+ Cot 



(3.3) 



Since f > j /, this yields the lower bound 

Ho 



m > 



Ci 
l + Cot 



(3.4) 
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Step 2: Longtime convergence in L^(il) x L^(0) x L^(ri). We now show the long- 



time convergence in Lp'. Estimates (2.8) and (2.15) respectively rewrite, for homogeneous 
boundary conditions, 



1 d^ 

2dt 



Gp\\u{t,-)\\l. + \\\T{t,-)\\l,)+G^ 



du 
dy 



(t,-) 



+ 



L2 



fr it 



L2 



(3.5) 



and 

1 d 

2dt 



dU 
dy 



(t,-) 



L2 2p 



d'^u 



dy^ 



it,-) 



L2 



<C2(\\f{t,-)h^\\{Vfr){t,-)\\l^ + \\u{t,-)\\l2 



where U is defined by (2.13). The evolution equation on r writes 



Af + /> = -/(r-.). 



(3.6) 



(3.7) 



We introduce the positive scalar e, to be fixed later on. We use the Cauchy-Schwartz and 
Young inequalities 



fir -f)\< \/f\f\ II ( V/(r - f) ) (i, 011^^ < ef\ff + ^ 



fir-f) (t,.) 



L2 



SO that, multiplying evolution equation (3.7) by r, we obtain 
Ad|fi2 



^^^it) + il-e)f\f\\t)<l-Jf(t,.)\\^. 



r-f)(t,. 



lL2 



The evolution equation on r — r reads 



A|(r-.) + ^(.-.) = -(/.-/.) + G|. 



(3.8) 



(3.9) 



Multiplying evolution equation (3.9) by r — r, integrating over il and repeatedly using the 
Young inequality, we find 



2dt 



(A||(r-r)(i,.)||iO + ^ll(r-r)(t,.)|lI. <cf||(/r)(t,. 



Il2 



+ 



dU 
dy 



itr) 



L^. 



SO that, using the uniform in time L°°-bound on /, 

^'^:(A||(T-f)(t,-)|li2) + ^||(r-f)(t,.)||i. 



2dt 



2i] 



<C3[\\fit,■)\\L^\\{^/fr)it,■) 



L2 



+ 



dU 
dy 



itr 



^2, 



(3.10) 



We introduce some positive scalars mi , ?tt,2 , m^ and the energy function 



E{t) =mi{Gp \\u{t, ■)\\l, + A ||r(t, ■)\\l,) + m2 
+ m3A||(r-f)(t,-)|li2+A|f(t)p, 



dU 
dy 



it,-] 



L2 



(3.11) 
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which therefore satisfy, combining (3.5), (3.6), (3.8) and (3.10), 



~{t) + {mi - C,m2 ||/(i, OIIli - C^m, ||/(t, ■)\\^^) || (//r) (t, 



+ {miGrj - C2m2Cp) 



du 
dy 



(t,-) 



V 
+ I 7r'^2 - Czm:iCp 



2 
L2 



9y2 



(i,-) 



L2 



+ I ^"^3 - ^ \\f{t, OIL^) ||(t - r) (t, Olli^ + (1 - 6)/|f |2(t) < 0, (3.12) 

where Cp is the Poincare constant. 

The coefficients mi,m2,m^ are chosen sufficiently large so that, for all time t > 0, 



every term in the left-hand side of (3.12) is positive. The conditions 
G 1 

— 1713 > — SUp||/(t, OIIloo, 



2p 



1712 > C^Tn^Cp, 

(G m G 
"^ J' ^ ,C2m2 sup \\f{t, Ollii + GsTUs sup \\f{t, •)||^„ 



are sufficient. Using in addition the lower bound (3.3) and the Poincare inequality, (3.12) 
becomes 

2 



IdE 



{t) + Co{\\u{t,-)\\h + \\{r-f){t,-)\\h + 



dU 
dy 



it,-) 



2 dt 

Using the triangle inequality 

mi I|t((, Ollia < 1 1f I" (<) + e ||(t - f ) (t,.,||^, , 



,.,+<^-)rTWi^i^<"^»- 



2mi 
we find, for t sufficiently large 



Ci 



IdE 1 ( e 

— ; — i — mm 1 — 2e, — , ^ 

2 dt A V 2miJ l + Got 



E<0. 



(3.13) 



(3.14) 



We take e < - and apply the Gronwall Lemma to (3.14) to obtain that E goes to zero in 
the longtime limit. In particular, we have 



hm l\\u{t,-)\\i2 + Mt,-)\\i2 + 



dU 
dy 



(t,-) 



i^ 



(3.15) 



Step 3: Longtime convergence in H^{Q) x L°°(il) x L°°{Q,). Combining (2.15) 

dy' 



and (2.17), using the uniform in time L°°-bound on / and the Poincare inequality on 



the spatial average of which is zero, we obtain 



d^ 
dt 



dU 
dy 



it, 



L2 



+ A|r(t,y)|2 + 



d^u 



<Go[Mt,.)\\l, + \\T{t,.)\\l,+ 



dy"^ 
dU 



it. 



L2 



+ \rit,y)\' 



dy 



it,-) 



L^. 
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We apply the Gronwall Lemma and use the convergence (3.15) that is uniform in space 
to derive 



hm ||T(t, ■ 



t— >oo 



0. 



Using the convergence of ||r(t, 
sufficiently large, 



l^oo, the evolution equation (1.1c) on / implies, for t 
df 



dt 



< 



1 



f. 



(3.16) 



This yields the convergence of \\f{t, ■)\\ioo to zero in the longtime. 

Additionally, using the definition (2.13) and (3.15), — — converges to zero in L^($7). 

1 r H r Qy 

This ends the proof. () 

We now turn to making precise the rates of convergence to the steady-state. We 
introduce the non-negative scalar 

/3 = meas{yGf]|/o(y)>0}. (3.17) 

By assumption in this section, we have /3 > 0. The following result establishes the con- 



vergence rates in function of /3. In Section 5.1, we will check using numerical simulations 
that these rates are indeed sharp. 



Theorem 3.2 Supply system ( 1.1 ) with homogeneous boundary conditions and initial con- 



ditions that satisfy (2.1) and /o ^0. The solution {u,t, /), the existence and uniqueness of 



which have been established in Theorem 2.1, satisfies the following convergence estimates: 
for any arbitrarily small a > 0, there exists a constant Kq, independent from time and 
there exists a time t^, both depending on the domain, the initial data, the coefficients in 
the system and a, such that, for all t > tQ, 

.It 



ll^(t,-)||^i + ||r(t,-)llL2<^a(l + t)~^^'-"\ 

where j3 is defined by ( |3.17 ) and for all t > to and y ^ 0. , we have 



1 



<f{t,y)< 



1 



j^^+{l + a){t-to) ^ + (l-a)(t-to) 

In addition, there exists another constant k^, such that, for all t > tQ, 



W{t,-)\\m + \\{T 
du 
'd-y^^ 



-r)it,.) 
r]it,-) 



L2</ta(i + t)-^-^(^-"\ 



<K„(i+r'-'('-"\ 



(3.18) 
(3.19) 

(3.20) 
(3.21) 



L2 



where the function r is the spatial average of t, defined as in (2.7). 



Proof. The proof falls in four steps. We first consider the fiuidity, then derive first 
convergence rates for the velocity and the stress. A study of the auxiliary function defined 



by (2.13) next allows to conclude on the convergence estimates (3.20) and (3.21). 

We fix e an arbitrarily small positive scalar, actually equal to — , where a is the constant 

that appears in the statement of the Theorem. The constants k^ depend on e and have 
value that may vary from one instance to another, the actual value being irrelevant. 
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Step 1: Convergence rate for the fluidity. In view of Theorem 3.1 \\T{t, 



and \\f{t, ■)\\^ao vanish in the longtime. This imphes that there exists a time Iq, such that, 



for ah y S n and t > to, the evohition equation (1.1c) on / leads to 



{l + e)f{t,y)<^{t,y)<-{l-e)f{t,y). 



(3.22) 



For all y E fi such that fo{y) > 0, we have /(to) v) > 0, as shown in Step [T] of the proof of 
Theorem 2.1 The equation ( |3.22 ) becomes, for such y and t > to, 

^ <f{t,y)<—, — , ,/ ,,, , , • (3.23) 



f{to,y) 



+ (l + e)(t-to) 



f{to,y) 



+ {l-e){t-to) 



and hence (3.19 ). This inequality is also valid for all y such that /o(y) = 0, that is f{to,y) 
0, and therefore for all y £ ft. 



Step 2: First convergence rates for the velocity and the stress. We first 
make more precise the lower bound on /. As /o is continuous, there exists a closed set Q^ 
such that foi^e) > and measjrie} = /3(1 — e). As shown in Step [I] of the proof of 
Theorem 2.1, we also have f{to, ile) > 0. Furthermore, as f{to, •) is continuous, we obtain 
/(to, ^e) > /«e- The inequality ( 3.23| ) thus becomes, for all y G il^ and t > to, 

f{t,y)> ^ ^ 



l + tK,+t 



It follows from f > f that, for all t > to, 



m > /3 



1 



1 



l + eKe + t 



(3.24) 



We now use the energy E introduced in (3.11). As ||/(t, OIIloo vanishes in the longtime, 
the coefficients mi,m2,m-i can be chosen arbitrarily small in (3.12 ), independently from e. 



for t > to sufficiently large. We insert (3.24) in (3.12) so that, for sufficiently large t, 

2 \ 



~{t) + Col\\u{t,.)\\l, + \\{r-f){t, 



\l^ + 



dU 
dy 



(t, 



L2 



Using the triangle inequality (3.13), we obtain, for sufficiently large t, 
IdE 



/3 
, + — min 1 
2 dt X ' 

and therefore, using that 1 — 2e < 



2e, 



1 



1 



2mi J 1 + e Ke + t 



-E< 0. 



2mi 



as mi is arbitrarily small, 



Ids /? 1 

2 dt ^ X^ Ke + i 



E<Q. 



Applying the Gronwall Lemma to ( 3.25| ), we find, 

dU 



Ht,-)\\l. + \\T{t,-)\\l,+ 



dy 



{t,- 



<^^(l + t)-2f(l-4.)_ 



(3.25) 



(3.26) 



L2 



where we recall that k^ denotes various constants. We have obtained (3.18) 
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Step 3: Convergence rate for the auxiliary function U. We recall that the 



function U is defined by (2.13). We first prove that U is more regular than claimed 



in (2.16). We rewrite (2.14) 



dU r]d^U 
dt p dy"^ 



I fy G 

T- / ifr - fT)dx + —u. 
A?? Jo Ar? 



(3.27) 



dU 



We deduce that — — satisfies the heat equation with a right-hand side in Lf ((to, +00), L ) 
oy 

dU 
and initial condition ^^(^0; ■) £ H^(Q,) at time to (up to a possible modification on a set 
oy 

dU 
of times of measure zero). Therefore, we have -—- £ Hl^^{(^to,+oo),L'^), so that 



UeHUito,+^),H}^] 



We next differentiate (2.14) with respect to t, insert (1.1a) and find, 



d^U 1] d' 



dU 



Xp dy^ 



dt"^ p dy"^ \ dt 
where / is the function defined by 

1 fy fdfT djr 



lit,y) 



>^V Jo 



dt 



dt 



L 



dx. 



(3.28) 



(3.29) 



(3.30) 



We now regularize I as follows. We consider a sequence of functions Im such that for all 
m, Im is infinitely differentiable from (to, +00) to L^(J7) and as m — )• 00, 



/„^ ^ / in Lf^ciito, +00), L^) 
Consider a solution Um G C°°((to, +00), H'^ D H^) to 

Gd'^Um 



d^Um V d^ fdUrr. 

dt"^ p dy"^ \ dt 



Xp dy"^ 



^m.' 



(3.31) 



(3.32) 



Equation (3.32) has been studied in [6, 7J. Inspired by arguments from these references, 
we introduce the energy functions Hm and Fm depending on a constant 5 G (0, 1) to be 
determined later 



Hm{t) 



dUrr 



dt 



it,' 



L2 Ap 



dU„ 



dy 



it, 



+ ^5 



L2 



dUrr 



dy 



it,-) 



+ S 
L2 Jn 



dU„ 



and 



Fm{t) 



V ,2G 



Xp 



+ 2'^6^ 



d dUrr 

dt dy 



d dUrr 

dt dy 



it, 



G j-2 



dt 

d fdUra\ dU, 



,.^'fp' J^KdtKdyJ dy 



Um ) {t, ■) 



it,-) 



2 

-5 

L2 


dUm 
dt ^ ' ^ 


2 

L2 




^ ( 

+ 5- 




dy ^' 


,•) 


2 

L2 



df^ 



it. 



2 

L2 
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We multiply (3.32) by 



dUrr 



+ 5Um + 5' 






, integrate over and find 



We use the Poincare inequality and choose 6 sufficiently small, depending on the domain 



and the coefficients in (1.1) such that, for suitable constants ci, C2 and C3, 
' d dU„ " 



c\ 



and 



dt dy 



'it,^ 



+ 



L2 



dU„ 



dy 



it,-) 



L^. 



< Hm{t) < C2 



d dU„ 



dt dy 



it, 



+ 



L2 



dU„ 



dy 



it,- 



L^. 



Fm{t) > C3 



d dU„ 



dt dy 



it,-) 



+ 



L2 



dUr, 



dy 



it,-) 



+ 



L2 



d'^Un 



df^ 



it,-] 



(3.33) 



(3.34) 



Using the upper bound in (3.33) and (3.34) and the Young and the Poincare inequalities, 
we obtain 



diiYi 
~dt 



it) + CH^it)<\\Im{t,-)\\L^- 



We multiply the above equation by e , integrate from to to t and find 

ft 



Hm{t)e^' < Hm{to)e^'' + / ||/™(s, Olll^ e^^ds 

Jto 



(3.35) 



(3.36) 



Equation (3.29) is linear so that by (3.28) and (3.31), we can pass to the limit m, — t- 00 



in (3.36) and find, for all t > to, 

H{t)e^' < F(to)e^*o + / \\I{s, ■)\\l, e^'ds. 
where H is defined as Hrry with U instead of Um ■ 



(3.37) 



The study of (|3.29|) reduces to the understanding of (|3.37|). We now make precise 

dfT 



the behaviour of / or more precisely at the one of 



and (1.1c) to ffiid 



dt 



We combine equations (1.1b) 



f^/OnD-^/^-z^-^M/^) 



(3.38) 



Multiplying the evolution equation (3.38) by /r and integrating over Q yields 



1 

2 


dt ^'' ^ 


L2 Jn\^ 

( 

<coii/(t,-)llioo 
v 



i-./ + fM)/V + f(A| 



Q 



du 
dy 



(t, 



L2 



+ \\T{t 



,-)\\l^ , 



(3.39) 
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where we have used the L°°-bounds on both r and / and the Cauchy-Schwarz inequaUty 
to derive the second Hne. Inserting (3.23) which gives the convergence in - of \\f{t, ■)\\ioo 
and (3.26), equation (3.39) implies 



dfr, 



dt 



t,-] 



<K,(l + t)-2-2f(l-4^). 



(3.40) 



L2 



dfr 



Since the L -norm of —- — controls the L -norm of /, we insert (3.40) in (3.37) so that, 

dt ^ ' ' ' 

for all t > to, 



H{t)e^^ < H{to)e^^° + k^ f — 

J to (1 



,Cs 



2+2f(l-4e) 



-ds. 



(l + s) 
Moreover, for q > 0, for all t > tQ, we integrate by parts to obtain 



t ^Cs 



-ds < 



to 



t gCs ^Ct 

ds 



{l + s)i - C{l + to)Jt(l + s)'i C{l + t) 



We insert ( |3.42| ) with g = 2 -h 2^(1 - 4e) in ( |3.41[ ), so that for t sufficiently large 

H{t) < K, 



1 



(l+t)2+2f(l-4.)' 



Using the lower bound in (3.33), we have therefore obtained 

2 



dU 
dy 



(t, 



L2 



<K,{i+tr^-^-^^'-'^\ 



(3.41) 



(3.42) 



(3.43) 



(3.44) 



Step 4: Convergence rates (3.20) and (3.21). Using (3.10) rewritten as 

~{M\ir-f){t,-)\\h) + ^\\{T-f){t,-)\\l, 



<C 



dU , 
dy 



t,-] 



+ ll/(i,-)llio.||r(t,.)||i. , 



and convergence estimates (3.23), (3.26), (3.44), we obtain 



{r-f){t,.)\\l,<.,{l + tr'~''^('-'^\ 



and eventually 



du 
dy 



it,' 



L2 



<K,(l + t)-2-2!(l-4e)^ 



We thus obtain (|3.20|) with a = 4e and conclude establishing (|3.21|) as follows: we return 

dfr 



to (3.39) and improve the convergence estimate for — - — , namely 
^ ' ot 



dfr 
dt 



it,-) 



L2 



</.,(l + t)-4-2^(l-^^). 
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This implies, mimicking (3.41) and using (3.42) with q = A + 2^(1 — 4e), that for t 
sufficiently large, 



dU 
dy 



it, 



L2 



<Keil + t) 



_4-2f(l-4E) 



that is (3.21) with a = 4e. 



(3.45) 





3.2 Case fo = 

In the case /o = 0, / vanishes for all time. System (1.1) then reads 

du 



dt 
dr 
dt 



d'^u dr 



dy' 



G 



du 
dy' 



(3.46a) 
(3.46b) 



The existence and uniqueness of a regular solution to (3.46) is easy to establish. The 



longtime behaviour of system (3.46) is now made precise. 



Theorem 3.3 Supply system (3.46) with homogeneous boundary conditions. Consider a 
solution (n, r) in the space 

(C7([0, +<x)-H^) n LL([0, +oo); H^)) x C([0, +oo); H^) 

Then, the solution converges exponentially fast to the steady state (0,7^) in H^{Q) x L^(r2) 
in the longtime: there exist two constants C , independent from time and initial data, 
and Co, independent from time, such that, for t sufficiently large , 

^''' - +||r(i,-)-T^IL2<Coe-^*. (3.47) 



dy 



it,' 



L2 



Proof. We perform the same manipulations as those used to obtain equation (3.29) in 



Step [3] of the proof of Theorem 3.2 Since we deal here with the case / = 0, we have 1 = 
in ( |3.29[ ). We have proven that studying the longtime behaviour to (3.29) amounts to 
proving (3.37). We therefore find, for t sufficiently large, 

dU , 



dy 



itr] 



<Coe 



-ct 



L2 



We next differentiate equation (3.46a) with respect to t and insert (3.46b) to obtain 

d'^u T] d'^ / du\ G d'^u 
dt"^ p dv"^ \ dt 



0. 



pX dy"^ 
The function u satisfies the same equation as U and thus has the same convergence rate. 



Applying the Gronwall Lemma to (3.9) therefore implies, for t sufficiently large, 



Integrating (3.46b) over Q,, we have 



r) (t, ■ 



X — T 

dt 



lL2 



0, 



(3.48) 



so that r(t) = tq for all times. We thus have the convergence estimate (3.47) 
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4 Longtime behaviour for non-homogeneous boundary con- 
ditions in a simple case 

In this section, we study the longtime behaviour of the system ( |1.1[ ) supphed with non- 
homogeneous boundary conditions n(t, 0) = and u(t, 1) = a (where a is a constant scalar 
different from zero and chosen positive, without loss of generality, a > 0). 



We denote (uocToo, /oo) a stationary state to the system (1.1). We only consider the 
simplified case 



/oo > everywhere. 



(4.1) 



The only stationary state that satisfies (4.1) is made explicit in subsection 4.1 In 



subsection 4.2, we show convergence in the longtime to this stationary state for small 



initial perturbations. In subsection 4.3, we study the longtime behaviour for initial data 
that satisfy /o > without any smallness conditions, but only in a simplified case that 



reduces system (1.1) to a system of ordinary differential equations. 



We do not state any result for the convergence to stationary states when fluidity 
vanishes on some part of 0. 



4.1 Stationary state 



The following lemma makes precise the stationary state that satisfies the condition (4.1 ) 



Lemma 4.1 (Stationary state) Supply system (1.1) with non-homogeneous boundary 
conditions Uoo{0) = and Uoo{l) = a > 0. The unique stationary solution {uoo,Too, foo) 
in 



{^H^{Q,)) satisfying (4.1) reads 



yuoc 5 Too 1 /oo )\y) 



ay, 



VI + 4.uCGa + 1 VI + 4i^CGa - 1 



2e 



2u 



(4.2) 



Remark 4.1 It is easy to extend the above result to stationary solutions (uoo, Too , foo) 
in H^{Q) X L°°(f2) x L°°{Q) that satisfy /oo ^0. Introducing floo = {v ^ ^,/oo(y) > 0} 
and /3oo = meas(Ooo), the set of such stationary solutions reads 



dupc 
dy 



.TooJoo (y) 



a -l + ^TL 

Poo J^ , 



on J7oo 

on i7\r2oo, 



1 



with Ti = -- ( 1 + ^1 + Aiy^Ga/Po 
2^ \ 

Proof. The stationary states (uoo, Too, /oo) : ^ 



of the system (1.1) that sat 



isfy (|4.1|) are solutions of the following system 

= rj 



d'^u, 



oo "-"oo 



= G 



dy"^ 
dy 

i + Ckc 



dy 

joqToo, 



(4.3a) 
(4.3b) 
(4.3c) 
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We now show that such a steady state is unique and exphcitly identify it. Since Tqo € 



H^{^), (4.3a) shows that Uoo belongs to H'^{^). We integrate (4.3a) and (4.3b) over f] 
and obtain 



K = 1]-- hToo, 

oy 
where if is a constant and, using the boundary conditions on Uoo, 



We combine (4.3b) and (4.4) to obtain 

(^/oo + l) Too = if 



(4.4) 



(4.5) 



(4.6) 



so that Too has the constant sign of K. Equation (4.5) then implies that Tqo, thus K are 
positive. 



We now claim that Tqo is constant over $7: inserting (4.3c) in (4.6), we obtain that r^ 
satisfies 



Too 1 + ??- 



-l+^To 

Gu 



K. 



It is easy to see that this equation has a unique positive solution Too. It follows from (4.4) 

UVL 

that ^i; — is constant throughout $7 so that, using the boundary conditions, Uoo{y) = ay. 



dy 



We rewrite equation (4.3b) as 



to find the value of 



Ga = Try 



Too = 7^(1 + v'l + 4j^eGa). 



_1 _)_ ^T^ 

The stationary state reads ( ay,Too, — ) , that is (14. 2|) 







4.2 Longtime behaviour with smallness assumption 



The following theorem states the convergence in the longtime to the stationary state (4.2) 
for small initial perturbations. 



Theorem 4.1 Supply system (1.1) with non-homogeneous boundary conditions u{t,0) = 



and u{t, 1) = a > 0. Consider the solution {u, r, /) the existence and uniqueness of which 



have been established in Theorem 2.1 and the associated stationary state (uoo,Too, /oo) 



defined by (4.2). There exists e > (sufficiently small so that at least tq and /o are 



positive ), such that, if the initial data (uq,tq, fo) for (1.1) satisfy 



\uo — UooWh'^ + Iko ~ TooWloo + ll/o — fooWioo < 6 , 
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then the solution {u,T,f) of system. (1.1) converges, as t goes to infinity, to {uoo,Too, foo) 
inH^{Q) xL°°{n) x L'^{n). 

More precisely, there exist a constant C independent from e, time and initial data and 
a constant k^ independent from time such that, for t sufficiently large, 



\u{t,-) -Uoo{- 



\H' 



+ \\r{t, 



+ ||/(t,-)-/oo(-) 



< KfC 



<C~e)t 



(4.7) 



Remark 4.2 It is indeed possible, under the same assumptions, to prove that both r and f 
converge to zero in H^(Q.) and not only in L°°(il). The proof is more tedious. We omit 
it here for brevity and refer to 11]. 



Before we get to the proof , we note that we will return to system (1.1), and not (2.4) 



since of course boundary conditions will play a crucial role throughout the section. We 



also rewrite system (1.1) as 
du 



dt 
dr 

m 

dl 
dt 



d'^ d 

'n-^{{u - Uoo) + Woo) + fl"(('^ ~ '^°°) + "^oo)' 

d 

G^-{{u - Uoo) + Uoo) - ((/ - /oo) + /oo)((t - Too) + Too), 

dy 



(-1 + el(T - Too) + Too|)((/ - /oo) + /o 



K(/-/oo)+/o 



To lighten the notation, we henceforth denote {u, r, /) instead of (n — Uoo,t — Too, f — fi 
and consider 

du d'^u dr 



P 



dt 



9r _ du 
dt dy 



dy"^ dy ' 

JooT Too J JT, 



(4.8a) 

(4.8b) 
uy 

§[ = -^{f + /oo)'/ + af + 2/oo)/r + UIt, (4.8c) 

supplied with homogeneous boundary conditions on u and initial data that satisfy 

II l|2 I II ||2 I \\ p \\2 ^ 2 

lPo||/fi + IrollL^o + II/oIIloo S e • 
Proof. The proof is divided into three steps. The first step establishes a priori 



estimates on system (4.8). In the second step, we show that the solution remains small 



for sufficiently small perturbations. In Step[3j we show that, still for small perturbations, 
the solution converges to the steady state and that the rate of convergence is exponential. 
As in the previous proofs, C and k^ denote various constants the value of which may 
vary from one instance to another, the actual value being irrelevant. 

Step 1: A priori energy estimates. We argue as in Step [2] of the proof of Theo- 



rem 



2.1 We multiply (4.8a), (4.8b) and (4.8c) respectively by ti, r and /, integrate over fi 
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and find 



P d ,, , , ,,2 



'^ ^ II U M|2 



(/ + /»)T=(«, ■)^a j At.-)-T^ I /t((, .), 

\2 j:2i 



~\\f{tr)\\l^+y JU + fooYf\t,.) = U'L jjr{tr)+i jj^f + 2f^)fMtr)- 



Combining these estimates leads to 



ld_ 
2dt 



(^pG \\u(t, Olli. + A \\r{t, Olli. + ^ \\f{t, Olli.) + r?G 



dy 



(t,-) 



L2 



+ /oo||r(i,-)|li2 



+ ^ii/(t,-)iii.-c 



e 



L^{L°°) 



.(L^) + 2/oo)(iiT(t,-)iii2 + ii/(t,-)iii2)<o, 

(4.9) 



using the L°°-estimate on [0, T] on / established in the proof of Theorem 2.1 and the 
Young inequality. We now use Step [3] of the proof of Theorem |2.1| and more precisely the 
estimate (2.15) on U defined by ( |2.13| ). We have 



dt 



dU 
dy 



it,- 



+ - 

L2 P 



d^u 



dy^ 



it,- 



L2 



<C \\u(L-) 



Il2 



+ 



2 

L^(L°°) 



\T{t 



,-)\\h 



(4.10) 



Now that we have estimates in Sobolev spaces, we turn to point wise estimates on r 
and /. We refine our argument in Step [4] of the proof of Theorem 2.1 We rewrite the 
evolution equation (4.8b) as 



,dT / ^ G\ BU G 

A^+/ + /oo + -r = G—- + -T- T^f, 
at \ rj J oy r] 



multiply it by r, apply the Young inequality and obtain 



^^^H^^^-^l'i^i'^" 



dU 

dy 



it,' 



+ \\Tit,-)\\l2 -T^fr. 



Similarly, we multiply (4.8c) by / and find 
1 d 



2dt 



i/r + K/ + fooYifr = af + 2/oo)r|/i^ + e/4/r. 



dU 



(4.11) 



(4.12) 



We combine (4.11) and (4.12) and use the Poincare inequality on — — , the spatial average 

dy 



of which is zero, to obtain 



A d 2 , •^oo ^ I f |2 , / f , / I "^ \ U|2 



2dt' ' 2(fldt' 



G 

2ri 



+ JW ( K/ + fooY - mfWL^iL-) + 2/oo) ||r||^. 



< 






'{L°°) 



l/p 



1] 



dy"^ 



it,-) 



L2 



+ Mt, 



Il2 



(4.13) 
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Step 2: Smallness of the solution for small perturbations. We now prove 
that, for e € (0, 1) to be fixed later on and an initial condition satisfying 



I l|2 ill ||2 ill^l|2 ^2 



we have, for all time t > 0, 



iKi,-)ii^i + iiT(t,-)iiioo + ii/(t,-)iiio.<6. 

We argue by contradiction and suppose 

Tm = inf {t G M+ I (\\u{t, Oll^i + ||r(t, ■)\\l^ + ||/(i, ■)\\l_ 



(4.14) 



(4.15) 



> e > is finite. 



For all t < Tm, we use the estimates from the previous step. For e sufficiently small such 



we have, integrating (4.9) from to t, 



that all the terms in the left-hand side of (4.9) are positive (this gives one condition on e 



pG\\u{t,.)\\l.+X\\r{t,.)\\l. + ^\\f{t,-)\\l2 



^ /o' V° 



du , 



L2 



I J'^ II / M 1 2 I ''"oo ^ 



\\f{sr)\\lAds<Ce\ 



Integrating (4.10) from to i then yields 

dU 



dy 



(t, 



+ 

L2 JO 



d^u 



dy^ 



i^r) 



ds < Ce\ 



L2 



We now integrate (4.13) from to t and get 



Try 



— \ |2 _|_ '°° I -fl2 

2 1^1 ^2^1 



\f\'<Ce'. 



(4.16) 



(4.17) 



For all t<TM, ( |4.16D and ( |4.17D imply 

du 



dy 



it, 



L2 



+ Mt,-)\\U + \\f{t 



' 7IIL°° 



<Ce' 



Choosing e sufficiently small such that Ce^ < e (which gives another condition on e) 



contradicts the definition of Tm, and so Tm = oo. It follows that (4.15) holds for all 
time t > 0, the solution remains small. 



Step 3: Convergence to the stationary state. We now prove that, if the initial 



data satisfy (4.14), then the solution converges exponentially fast to the stationary state 



in the longtime. For t sufficiently large, (4.9) implies that 



\n{t,-)\\l, + \\T{t,-)\\l, + mt,.)\\l,<K,e-'^^-^^K 



(4.18) 
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Adding ( |4.10| ) multiplied by 2p to ( |4.13[ ) leads to 

2 






2p 



at/ 

9y 



it,-) 



L2 



.^lp.^i/pU, 



52 c/ 



ay2 



it,-) 



+ if + foe + ^)\r\' 

L2 ^^ 






+ 7^- K/ + fooY - mfh^iL^) + 2/00) l|T||ioo(^^) I/I 



< 



c^ (h(t, oiii^ + iiT(t, oiii. + ii/(t, oiii^ + ii/iii5c(^..) iiT(t, oiii^) 



dU 



We use the Poincare inequality on — — , the spatial average of which is zero, apply the 

oy 



Gronwall Lemma, insert (|4.18|) and find 
dU 



dy 



it,-) 



+ l|r(t,-)|lio. + ||/(t,-)llio.<K.e-(^-^)*. 



L2 



This convergence estimate is equivalent to (4.7) and we have exponential convergence. 



♦ 



4.3 Longtime behaviour without smallness assumption (simplified case) 



We now examine the longtime behaviour of system (1.1) supplied with not necessarily 



small initial data iu(),To, fo). We are unable to prove a general result and focus our 
attention to the particular case where the initial condition is uq = ay (a positive constant). 
To = constant = tq, fo = constant = /o > 0. In such a case, a substantial simplification 



occurs. Indeed, (|1.1|) reduces to the following system of ordinary differential equations: 

(4.19a) 



( dr 



dl 

[ dt 



-JT + Ga 

i-l + i\r\)f-yf. 



(4.19b) 



supplied with initial conditions tq, /o G M with /o > 0. 



System (4.19) has a unique steady state such that /oo > and it reads 

(,'''00 1 /oo j - 



yjl + Av^Ga + 1 VI + AviGa - 1 



2C 



2v 



(4.20) 



Indeed, such a steady state (too,/oo) satisfies /ooToo = Ga (so that Too > 0) and z//oo = 
— 1 + ^Too. Combining these equations implies (— 1 + ^Tqo) Tqo = vGa. This equation has 



a unique solution given in (4.20) 



In addition, we introduce the condition 



^-^«--)M(^-i 



<0, 



with 



3Ga 



a = mm 



\ Gau + 4roo ' 



VI + Ay^Ga - 1 
3z^ 



(4.21) 



(4.22) 
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We are unable to perform our proof without this additional assumption. The numerical 
simulations we perform (see Figure p| however show convergence holds even when (4.21) 
is not satisfied. 



Theorem 4.2 Supply system (4.19) with initial conditions tq, /q G M with /q > 0. Then 



the solution (r, /) remains bounded. 



In addition, under assumption (4.21), the solution {t, f) converges to (roo,/oo) in the 
longtime and the rate of convergence is exponential: for t sufficiently large, 



Ht)- Tool + \f{t)- fool <Coe-^^\ 
where Cq is a constant independent from time and Cr reads 

(^ + ^fl)-lVA, if A>0, 

with 



Ot. 



if A < 0, 



(4.23) 



(4.24) 




+ ^/o 



4 I jfoo + T-Too 



(4.25) 



Proof. The proof is divided into seven steps. Step [T] introduces simplifications on 
the initial data and the system, that are not restrictive for the longtime behaviour. Some 
notation is given in Step [2| A lower bound on / is derived in Step [3] and is used in Step |4] 
to prove that the solution is bounded. Further restrictions are made in Step [5] still without 
loss of generality. Step [6] establishes the convergence, which is proven to be exponential 
in Step[7[ 

We consider until Step |4] the maximal solution to (4.19) although the solution a pos- 
teriori exists for all times because of boundedness. 



Step 1: Simplifications on the initial data. We show that r and / solution 



to (4.19) remain positive, possibly after some time for r. We first remark that, since /o > 



0, / > for all times, arguing as in Step[T]of the proof of Theorem 2.1 On the other hand, 
if T < on some time interval, evolution equation (4.19a) thus implies that r increases 



strictly on this time interval (recalling that a > 0). Hence, there exist a time Tq such 

that t{Tq) > 0. Moreover, for all t > Tq, t remains positive (since if r is zero at one 

(It 
time Ti > Tq, -^(Ti) = Ga > 0, which is in contradiction with r > for t < Ti.) 

For the purpose of studying the longtime limit, we may always consider, without loss 
of generality, the system 

dr 



A— = -/r + Ga 

-l + ^r)/' 
supplied with positive initial conditions tq, /o- 



dt 
dl 

dt 



^f 



(4.26a) 
(4.26b) 
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Step 2: Some notation. We consider the tiiree subdoinains: 

^i = |(r,/)|/>^^^, f<a\, 



A2 = UT,f)\f<^-^^, f<a}, 



A3 = {{T,f)\f>a}, 
where we recall that a is defined by ( |4.22[ ). We also introduce their intersections: 

ri3 = |(T,a)|r<^^^^ 



ri2 = <!(T,/)|/ = ^^,/<a 



r23 = Ut,<7)\t> 

See Figure [T] for a graphical description. 




Figure 1: Notation on (0, +00) x (0, +00) 



Step 3: Lower bound on /. We now establish a lower bound for the fluidity / in 
each domain. In the cases (tq, /o) G ^2 or (tq, /o) € ^3, we have 

/ > min{/o,cr}. 



The case (to,/o) G Ai requires more developments. The evolution equations (4.26a) 



and (4.26b) respectively rewrite 



d 1 



A d 

2dt l^~ I 



1 - ^T + uf, 
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We combine these two equations and obtain 
^d (Gal If 4^^' 



= -3Ga + Gavf - ^fr^ + 4/r 

< -3Ga + {Gau + 4roo) a 
<0, 



(4.27) 



where we have used firstly that < / < u and r < Too in Ai and secondly ( 4.22[ ). 
Integrating (4.27) yields 



Ga I Ga fill I f i^a + I 4 



ua + 1 
using that r < — - — . We introduce 



m 



f 



111 A^ fiya + l 4 
max <(—,-> + — — 7 h 



/oVJ 2Ga V e e 



(4.28) 



which is therefore a lower bound for / in the region Ai. This lower bound also holds for 
initial conditions that belong to A2 and A3 and we thus have, for all t > 0, 



fit) > ruf, 



(4.29) 



with ruf defined by (4.28) 



Step 4: Boundedness. The purpose of this step is to prove that the solution (r, /) 
remains bounded. 



Applying the Duhamel formula on (|4.26a|) yields 

Ga /■* _ r^- lldlri.' 
'0 



r(t) = e" Jo ^'^Vo + ^r / ^~ ^ ^^' 

A Jo 



so that, using the lower bound (4.29) on / 



r(t) < e A To + 






Therefore, r is bounded, and there exists a time to such that, for all t > to, 

Ga + l 



r{t)< 



m 



f 



(4.30) 



We now turn to the boundedness of /. We introduce M^- = and Mf = — (— 1 + ^Mr). 

We will show that, for all t > to, 



/(t)<max(/(to),M/) 



(4.31) 
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df 

We distinguish between two cases. Let us first suppose that ^-(^o) ^ 0. In this case, 

/(io) < -(-1 + ^T{to)) < Mf because of ( |4.30[ ). Moreover, for all t > to, f{t) < Mj. 
Indeed, by contradiction, if 



h = inf {t > to, f{ti) = Mf} < +00, 



df 
then, by continuity, -^(^i) ^ 0. On the other hand, we have 



dl 
dt 



(ii) = /'(ti)(-l + er(ti)-z^/(ti)) 



<0, 



uMf 



hence the contradiction. 

df 
In the other case ^(^o) < 0, / strictly decreases until (possibly) equality occurs at a later 

time ^3 I ^-(^3) = ) ) which leads to the previous case with t^ instead of to- In any case. 



we have obtained (4.31) for all t > to- 



Step 5: Further simplifications on the initial data. Table 1 first summarizes 
how (r, /) behaves when it touches an intersection line. We use Table 1 to show that the 



starting line 


t>f 
at 


at 


entering region 


ri3 

ri2 
r23 



+ 


+ 
+ 


A2 
^3 



Table 1: Motion on intersection lines 



solution enters region A^ at some time. 

ar 

In region Ai, we have — - > 0, so that there does not exist any periodic orbit inside 

at 

region Ai. There is also no steady state in this region. Using the Poincare-Bendixson 



Theorem on the bounded solution of ordinary differential equation system (4.26), the 



solution leaves region Ai at some time. According to Table 1, it enters region A2. 

df 
Applying similar arguments on region A2 where — - > 0, the solution enters region A3 at 

at 

some time. 

We can therefore restrict the studying of the longtime limit to initial data (to,/o) that 

belongs to region A^, , without loss of generality. 



The bounds (4.29) and (4.30) become 







(4.32) 
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and 



-("^(""^^'U^i; 



ua 



1 4 



(4.33) 



Step 6: Convergence. We introduce 

G(t,/) = - 



> + i«" 



and 

We have 

dG{Tj) , 5F(t,/) 



FiT,f) = i-l + Cr)f-uf. 



dr 



+ 



9/ 



</ 



A 

1 

A 



2 + 2^r - 3z^/ 



-^«>--)^^(^n 



using (4.33) and the positivity of /. Because of our assumption (4.21), the right-hand 



sides F and G of (4.26) satisfy 



dGirJ) , dF{Tj) 



dr 



df 



<0. 



According to the Dulac Criterion, there does not exist any periodic orbit for (4.26). Since 
it has only one steady state (tqc/oo), the solution converges to it: 



lim(|r(t)-Too| + |/(t)-/o 

I— ^-oo 



0. 



Step 7: Exponential convergence. Now that we have convergence to the steady 



state, we can use linear stability. System (4.26) linearized around the stationary state (tqo, /o. 
reads 



dt\ fl 



A A 



n 
fl 



The eigenvalues of the associated matrix depend on the sign of A defined by (4.25). 
If A < 0, the eigenvalues are complex and their real part is —\\\+ i^f^ ) • If A > 0, the 



eigenvalues are real negative, the smaller one in absolute value is — ^ ( a + ^foo ) + 2"^^ • 
The real part of the eigenvalues gives the rate of convergence and hence of values of Cr 



in (4.24). 



♦ 
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5 Numerical results 

In this section, we present numerical simulations that complement the theoretical results 
on the behaviour of the previous sections. 

We simulate numerically (1.1) in the interval il = [0, 1] and the interval [0,T] for T = 
10000. The system is supplied either with homogeneous boundary conditions or non- 
homogeneous boundary conditions u(t, 0) = and n(t, 1) = a for all time t G [0,T]. In 
the latter case, we take a = 1. As for the initial conditions, we take sinusoidal functions 
for all three fields. The values of uq oscillate between —0.002 and 0.002 for homogeneous 



boundary conditions and between and a otherwise, 
between —0.5 and 0.5. 

We use the following set of physical parameters, 
viscosity rj = 1 so that the Reynolds number is low. 



coefficients ^ and u are equal to one. 
otherwise stated. 



The values of tq and /o oscillate 

The density p = 0.001 and the 
The elastic modulus G and the 
The characteristic relaxation time A is 0.5 unless 



92 



[Un 
(/n 



Un-lj 



V- 



Tn-l) = G 



dUn 



+ 



d 



Tn-l 

dy 



System (1.1) is solved using a constant time step At 
scheme: 

At 

At 

1 

At 



Jn—lTn 



/„_l) = (-l + 



dy 

l + C|rn|)/n 



0.005 with the following time 

(5.1a) 
(5.1b) 
(5.1c) 



ifn - i^/n-l/, 



For the space variable, we use linear PI finite elements for u and piecewise constant finite 



elements for both r and /. Note that, in contrast to the approximating system (2.24) we 
used for our theoretical proof, we take t„_i instead of r„ in the right-hand sides of (5.1a) 
and (5.1b). This allows us to solve each equation separately. This choice is made for 



simplicity. Other approaches could have been employed. For our tests, we use elements of 
constant size h = 0.002 and perform the computations using ScilabfTT]. 

5.1 Homogeneous boundary conditions 

We first focus on the homogeneous boundary conditions on u considered in Section [3} The 
case /o = 0, that implies / = for all times, is uninteresting numerically. We therefore only 
show results for /o ^0. In this case, we have convergence to the stationary state (0, 0, 0) 
as proven in Theorem |3.1[ The convergence estimates are established in Theorem |3.2[ We 
recall the parameter 

^ = meas{y£n\fo{y) > 0} 
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and these convergence rates: for a arbitrarily small, there exist various constants Kq, and 
a time to ! such that, for all t > Iq, 



1 



+ (l + a)(t-to) 



\r{t,-)\\L2 <Ko,il + t)—^ 
<f{t,y)< — 



-f(l-a) 



/(to,?/) ^^"^"^^'^ '^u; f(^to,y) 

Mt,-)\\m + \\{T-T){t,-)h2<K^il+t) 



+ {l-a){t-to) 



du 



(t,-) 



L2 



< Ka{l + ty 



-1(1-") 



(1-a) 



Note that the last three estimates are exactly the same as in Theorem 3.2, the first 



estimate is an immediate consequence of (3.18) and (3.23). We now check that these 

1. The evolutions of 



estimates are sharp. We begin with the case /o > on il that is /3 
\\r{t, OIIl^ , \\f{t, OIL^ , Mt, Oll^i + ||(r - r) (t, -^^^ , || (^r/|^ + r - 



it,-: 



are repre- 



L2 



sented in Figure 2(a) We use a log-log representation. The slopes s, which correspond 
to a decrease as t^, are fitted on the numerical results and indicated on Figure 2(a)[ the 
numerical convergence rates, obtained with A = 0.5, are in good agreement with the esti- 
mates. 

We next consider cases where /o = on some part of the domain. In Figure 2(b)| we show 
simulations obtained with different values of /3. For each simulation, that is for each value 
of (3 considered, the convergence rates are fitted and represented as a function of f3. The 
numerical and theoretical convergence rates s agree. 
We have extended these results to the other values of A than A = 0.5 and other values of 



the parameters p,ri,G,^,iy to check that the convergence estimates of Theorem 3.2 depend 
only on A and /3 and are indeed sharp. 




(a) 







\u\}ji+\r-T 


Z.2 





"^^^--^ 


T 
f 

\T]dyU + r - T 


L2 

L2 □ 






^^^~~~~^ 


2 




^^^-^.^^^^ 


3 




^^~^-^^ ^^ ^^^ 


4 




~^-~~-^ 



(b) 



Figure 2: (a) Time evolution in log-log scale for homogeneous boundary conditions; the 
points are the simulated trajectories; the lines and the corresponding slopes s are fitted, 
(b) Fitted convergence rates s for /3 = 0, 0.01, 0.1, 0.6, 0.9, 0.99; the lines are the theoretical 
convergence rates function of (3. 
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5.2 Non-homogeneous boundary conditions 



The longtime behaviour for non-homogeneous boundary conditions has been studied in 
Section[4| We consider only stationary states {uoo,Too, /oo) that satisfy /oo > everywhere. 



We have shown that such a steady state (4.2) is unique. We established in Theorem 4.1 



that we have convergence to this steady state for small perturbations. To have convergence, 
we of course need to assume /o > everywhere. We observe numerically that no other 
condition, and specifically non assumption on the smallness of the data, is required. We 
consider the perturbations {u—Uoo, t—Too, f—foo) to equilibrium and show that they vanish 
in the longtime, see Figure [3J The evolution is plotted in semi-logarithmic scale. The 
convergences of the various norms ||r(t, ■)\\i2 , ||/(t, ■)\\l2 , \\u{t, ■)\\fji + ||(t — f) (t, OIIl^j 

— — \- T — t] (t, •) are indeed exponential. 



\ % s = -0.40174 


\t 
1/ 




\ *%= -0-36749 






; 


\ '': 




°a 


„;3;„„„„..„„....„„...„ 




\ -.N^^KXXXXXXXXXX.XXXXXXXXX^ 



10 20 30 40 50 



Figure 3: Time evolution of the perturbation to equilibrium in semi- logarithmic scale for 
non-homogeneous boundary conditions ; the points are the simulated trajectories and the 
line and the corresponding slope are fitted. 

In section |4j in order to establish a result without any smallness assumption, we have 
considered a particular initial data that reduces (1.1) to the ordinary differential equation 



system ( 4.19[ ). We have obtained convergence to the stationary state (4.20) and explicit 
formula for the rate of convergence. Numerically, we observe convergence even when the 



condition (4.21), which was assumed for the proof of Theorem 4.2 is not satisfied. The 



time evolution is shown in the space of (r, /) and the convergence are represented in 



Figure 4(a) We check that the convergence is exponential as observed numerically in 



the general case of ( |1.1[ ) (see Figure^. Moreover, we compute the convergence rate and 
compare it to the theoretical rate Cr defined by (4.24). The evolution of the perturbation 
function |r(t) — Tool + \f{t) — foo\ is plotted as a function of time in semi-logarithmic scale 
in Figure 4(b)[ The first case A = 0.5 correspond to the case when the eigenvalues of the 



associated linearized system are complex, the expected value of Cr is 0.8090; the other 
case A = 0.1 is when the eigenvalues are real negative, the expected value of Cr is 1.7895. 
The theoretical and numerical value agree. 



Acknowledgements. The authors are grateful to Frangois Lequeux (ESPCI Paris) for 
many stimulating and enlightening discussions on the rheology of complex fluids. 
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(a) 
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Figure 4: For the system of ordinary differential equations (4.19), time evolution (a) in the 
space (r, /) for A = 0.5; (b) of the perturbation to equilibrium in semi-logarithmic scale 
for A = 0.5, 0.1; the points are the simulated trajectories and the line and the corresponding 
slope are fitted. 
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